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We study the intricate relationships between the dynamical scaling properties of electron wave 
packets and the multifractality of the eigenstates in quantum systems. Numerical simulations for 
the Harper model and the Fibonacci chain indicate that the root mean square displacement displays 
the scaling behavior r{t) ~ t^ with /3 = D^, where D2 is the correlation dimension of the multifractal 
eigenstates. The equality can be generalized to d-dimensional systems as fi — /d, as long as the 



electron motion is ballistic in the effective D^-dimensional space. This equality should be replaced 
by /3 < D2 /d if the motion is non-ballistic, as supported by all known results. 
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Spatial extension of eigenstates influences electron dy- 
namics and transport properties of solids. It is well 
known that extended eigenstates in periodic lattices re- 
sult in ballistic motion of electrons while localized eigen- 
states in strongly disordered systems lead to zero mo- 
bility. As a consequence, one expects metallic phase 
and insulating phase, respectively. In between these two 
extremes stand the critical eigenstates, i.e., the multi- 
fractal eigenstates, which are neither extended nor expo- 
nentially localized. Multifractal eigenstates exist in var- 
ious quantum model systems, such as the Harper model 
of two-dimensional (2D) electrons in the presence of a 
perpendicular magnetic field tight-binding models of 
quasicrystals j^-j^ , and disordered systems at the metal- 
insulator transition (MIT) |]^-|l3||. In recent years, a great 
deal of effort has been devoted to the analysis of the mul- 
tifractal nature of the eigenstates. Effort has also been 
made to establish intrinsic connections between the mul- 
tifractality of the eigenstates and the transport properties 
in such systems. As an example, the unusual resistivity 
increase caused by decreasing temprature or improving 
structural quality of quasicrystals has been attributed to 
the multifractal nature of the electronic states [^. De- 
spite all these efforts, a fundamental level of understand- 
ing on precisely how the multifractality of the eigenstates 
affects the electron dynamics is still lacking. 

The characteristics of electron dynamics is reflected 
by the time evolution of the corresponding wave packet, 
whose spreading can be described by the root mean 
square displacement r{t). For ballistic spreading, r{t) ~ 
t, while in the absence of diffusion, r{t) ^ for large t. 
Numerical simulations jlj] for the Harper model and the 
Fibonacci quasiperiodic chain showed that r{t) ~ t^^ 
with < /? < 1. Such a scaling behavior has also 
been found in other systems with multifractal eigenstates 
1^ ill • Oil the other hand, the energy spectrum of 
the multifractal eigenstates can be measured by the mul- 
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The second moment, D2, 



tifractal dimensions 

has been shown to determine the decay of the average 
probability, C{t), of a wave packet remaining at its ini- 
tial position: C{t) ^ t'^^ It is naturally tempt- 
ing to link the spectral multifractal dimensions, Z?^, with 
the dynamical scaling exponent, (3 pO|-p^. Indeed, ear- 
lier such efforts have led to a variety of results, includ- 
ing [3 < D'^ p = D^_, ||,|3, and /3 > Dl^/ut 
p5[ , where D2 is the correlation dimension of the eigen- 
states. We note that all these results were derived with 
the assumption that the energy spectra are singularly 
continuous. Scaling behavior of r{t) may also exist in 
higher dimensional systems whose spectra are not singu- 
larly continuous ||l^,|l6|,|l9| , but it would be even more 
difflcult to derive a unique relation between /? and D^. 

In this Letter, we use a new approach to understand 
the scaling behavior of electron dynamics. Instead of 
focusing on the energy spectra, we study how the multi- 
fractality of the eigenstates in real space affects the time 
spreading of the wave packets. Our study leads to sev- 
eral fundamentally important flndings. First, for the ID 
Harper model and the Fibonacci chain, our numerical 
simulations show that (3 is directly given by the correla- 
tion dimension of the multifractal eigenstates: /? 



Second, by interpreting as the effective dimension 
of the multifractal space for electron motion, we gen- 
eralize the above relation to d-dimensional systems as 
(3 = D2 I d. This equality holds as long as the spreading 
of wave packets in the -space is ballistic in nature, as 
demonstrated by applying to exactly solvable 2D and 3D 
quasiperiodic systems. Furthermore, when the spreading 
in -space is non-ballistic, we have (3 < D2 /d, which 
is in complete agreement with all the known results for 
higher dimensional (d > 1) disordered systems |p|-|l3|. 

We start with a ID lattice of size L = Na, where 
N is the total number of sites and a the lattice con- 
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stant. The Hamiltonian has the tight-binding form 
H = X '^1 +Sm,n*m,n|m >< n|, where 

e„ is the site energy, tm,n the nearest-neighbor hop- 
ping elements. Eigen-energies and eigenstates are de- 
noted as Em and ipniEm), respectively. The multi- 
fractal analysis of an eigenstate is based on the stan- 
dard box-counting procedure: divide the system into 
a number of boxes of linear size wL (0 < w < 1), 
and then calculate the box probability in the kth box, 
Pk = J2nGk A Spatial Correlation function 

R{w,E) is defined as R{w,E) = X^fcPfe- general, 
i?(u', E) ^ ui^s (-E) ^ where D^^E) is called the correlation 
dimension of the eigenstate. Extended eigenstates and lo- 
calized eigenstates are characterized by D^(£') = 1 and 
(-^) = 0, respectively. For multifractal eigenstates, 
< (E') < 1. Considering the fact that eigenstates 
contributing to the time spreading of a wave packet may 
display different multifractal behaviors, we introduce a 
spectral-averaged spatial correlation function, R{'w) = 
Jj ^{'^1 Era), and define a spectral-averaged correla- 
tion dimension D2 if R{w) ~ w^^ holds. The root mean 
square displacement for a wave packet initially located 
at site no is given by r(i, no) = VSn I" " '^op|0«(t)p, 
where the wave function can be obtained by 

(t^n{t) = Y.m'^lo^Em)tpn{Em)ex.v{iEmt). For inhomoge- 
neous systems, r(t, no) varies with no- In this case, one 
averages r(i, no) over different initial wave packets in or- 
der to obtain the dynamics for the whole system. The 
scaling behavior of the system-averaged root mean square 
displacement is defined as r{t) ^ , with < /3 < 1. 
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FIG. 1. Multifractal properties of eigenstates for the 
Harper model of system size A'^ = 3000. The solid line shows 
R{w) ~ with = 0.5. The dotted and the dashed 

lines show R{w,E) ~ (^) for energies Ee and Ec with 
D'^{E^) = 0.38 and D^{Ec) = 0.64, respectively. 



Hamiltonian parameters in the Harper model U are 
tm.n = 1 and e„ = 2 cos(27r(Tn), where a = (\/5 — l)/2 
is the golden mean. We find that the correlation func- 
tion R[w, E) for various eigenstates displays a scahng 
behavior R{w,E) ~ w'^^^^\ where < D'^{E) <1 and 
D2{E) is different for different eigenstate. As examples, 
we show R{w, E) for a system of = 3000 in Fig. 1 at 
energies E^ = 2.59751453292 and E^ = 0.00000066284, 
which respectively correspond to the upper edge and the 
center of the energy spectrum. The corresponding scal- 
ing exponents are D'!^{Ec) = 0.38 and D'^{Ec) = 0.64. 
We also find that the spectral-averaged (over the whole 
spectrum of —2.6 < E < 2.6) correlation function R{w) 
displays a scaling behavior R{w) ^ w^^ with D2 = 0.5 
as shown in Fig. 1. In Fig. 2, R{w) is compared with r(t) 
obtained by averaging over 500 initial wave packets. We 
find r{t) ~ with P = 0.5, in agreement with a previous 
study Q. We therefore have [3 ^ = 0.5. 
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FIG. 2. Scaling behaviors of r{t) ~ r (thick lines) and 
R{w) ~ ui^2 (thin lines) for the whole spectrum (solid 
lines), the center sub-spectrum (dashed lines), and the edge 
sub-spectrum (long dashed lines) in the Harper model with 
P^Dt = 0.5. System size is iV = 3000. 

We can also study how eigenstates in different spectral 
regions affect the spreading of wave packets by consider- 
ing (j)n{t) = J2E,„en'^noiETn)tpn{Em) exp{iEmt), where 
il denotes the spectral region. The corresponding r{t) 
is compared with R{w) obtained by averaging the eigen- 
states within fi. Both the spectrum and the eigenstates 
of the Harper model have self-similar trifurcating struc- 
tures 0-0]. Therefore, we expect that the spreading of 
the wave packets contributed by eigenstates in different 
sub-spectra follow the same scaling law. In Fig. 2, we 
show R{w) and r(t) for the center sub-spectrum —0.2 < 
E < 0.2 and the edge sub-spectrum 1.8 < -E < 2.6. It is 
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clear that R{w) for different sub-spectra is the same be- 
cause of the same muhifractal behavior of the eigenstates 
in different sub-spectra. Furthermore, r(t) for different 
sub-spectra in Fig. 2 displays the same scaling law with 
/3 = = 0.5, as expected. This indicates that the mul- 
tifractality of the eigenstates indeed plays a crucial role 
in the spreading of the wave packets. We note that the 
downward shift of r(t) shown in Fig. 2 for the sub-spectra 
is due to the decrease in the number of eigenstates con- 
tributing to the spreading of the wave packets. 

We find the same relation (3 = D2 to hold for the Fi- 
bonacci quasiperiodic chain Q. In this model, i^.j = 1 
and Ci takes —u and u arranged in a Fibonacci sequence. 
The averaged correlation function has the scaling behav- 
ior R{w) ~ with < D2 < 1 depending on u. 
The system- averaged root mean square displacement r{t) 
used to be compared with R{w) was obtained by an aver- 
age of 500 wave packets. As examples, we show and 
r{t) for the whole spectrum in Fig. 3 around u — 1.65: 
(3 = = 0.62 for u = 1; /3 = = 0.5 for U = 1.65; 
and (3 = =0.43 for u = 2. These three cases cover 
the super-diffusive, diffusive, and sub-diffusive regimes of 
motion in real space. 
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FIG. 3. Scaling behaviors of r(t) ~ (dashed lines) and 
R(w) ~ (solid lines) in the Fibonacci chain with site 

energies it = 1, 1.65, and 2. The corresponding exponents are 
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0.62, 0.5, and 0.43, respectively. System size is 



The relation /? — D2 observed in different models can 
be generalized to higher dimensional systems with multi- 
fractal eigenstates. From its definition, R{w) is in fact the 
average probability for an electron in a box of linear size 
wL. Because R{w) ^ w^^ , the effective number of boxes 
of the system is N' = i-^)^^ , instead of in order 

to ensure the normalization N' x R{w) ^ 1. From the 
standard definition of the volume dimension Dy, i.e., the 



number of the measured boxesx{wL)^^ = constant, we 
have Dy ^ D^. Thus can be viewed as the effective 
volume dimension for electron motion in the multifractal 
space (hereafter called the D2 -spa.ee). If the motion (i.e., 
the spreading of the wave packet) is ballistic in D2-spa.ce, 
the volume, N', visited by the electron is N' ~ t^^ . Be- 
cause the width of the wave packet r(t) is still measured 
in the d-dimensional real space, we have r(t) ~ t^^^'^, 
which immediately leads to the result P = D2 /d. 

As tests of the above generalized relation, we consider 
the exactly solvable models in higher dimensional Fi- 
bonacci lattices [ p"5|jl9| ] and the Labyrinth lattices [^ . 
We find that the equality /3 — D2 /d holds for all those 
models. As an example, here we only show the case for 
3D Fibonacci lattice, which is constructed in a cubic lat- 
tice with constant hopping elements tij — 1 and sepa- 
rable site energy = e^a; -I- e^y -\- e^^, where e„, e^^, and 
eiz are the ID Fibonacci site energies. The energy spec- 
trum is given hy E = Ex-\-Ey-\-Ez and the corresponding 
eigenstates are ipi^E, x, y, z) = ilj.^{Ea:)x ipyiEy) x tpz(Ez), 
where Ej and ipj{Ej) {j = x, y, z) are the energy spectra 
and the eigenstates for the ID Fibonacci chain, respec- 
tively. As a consequence, these eigenstates are multifrac- 
tal , and their spreading dynamics is non-ballistic 
in real space: r{t) ^ t^ where < /? < 1 takes the 
same value as that for the ID Fibonacci chain. We note 
that all the relationships between the energy spectra and 
the dynamics ]2C|-[25|] cannot give correct prediction for 
the spreading of wave packets in these higher dimensional 
quasiperiodic systems. In our approach, for the higher di- 
mensional Fibonacci lattices, it is easy to show that R{w) 
has the scaling behavior with D2 = dD'^'' , where 02"° is 
the multifractal dimension of the eigenstates in the ID 
Fibonacci chain. Therefore, we have /? = D^/d. 

Electron motion in the effective -space can also be 
non-ballistic. In such cases, the above equality should 
be replaced by /9 < D^^jd, because no motion should 
be faster than ballistic. One special example is the dif- 
fusive motion (/3 = 0.5) in the metallic regime of 3D 
disordered systems with extended eigenstates {p^ — 3), 
where /3 < £'2/3. As shown below, this inequality is also 
satisfied by many other known examples involving non- 
ballistic motion in the Dl^-spaee. At the metal-insulator 
transition of 3D disordered systems (in the absence of 
magnetic field), the eigenstates are found to be multi- 
fractal with different values of D^: D2 /3 = 0.57 ± 0.1 
g, D^/S = 0.48 - 0.6 §, L»|/3 = 0.57 ± 0.07 [§, and 
D^/3 = 0.5 ± 0.07 For such systems in a mag- 

netic field, we have Df/3 = 0.57 ±0.07 @, D^/3 « 0.5 
|ll|, and D:^/3 = 0.43 ± 0.03 |l|]. On the other hand, 
the corresponding electron dynamics was found to have 
the same scaling exponent (3 = 0.33 jl^, showing that 
13 < D^/d for all these cases (with or without a mag- 
netic field). Existing numerical results for 2D disordered 
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systems satisfy this inequality as well. For the 2D dis- 
ordered systems in a strong magetic field, eigenstates at 
the centers of the Landau bands are niultifractal with 
D2 /2 = 0.76 ± 0.03 and the corresponding spreading has 
a scaling behavior with /3 = 0.5 At the MIT of the 
2D disordered systems with spin-orbit interactions, one 
finds D|/2 = 0.84 ± 0.03 and /? = 0.5 It is evident 
that /3 < D^/d. 

In conclusion, we have shown that there exists an in- 
trinsic relationship between the dynamic spreading of 
electron wave packets and the multifractality of the eigen- 
states in quantum systems. For the ID Harper model 
and the Fibonacci chain, we have (3 = D2, where /? is 
the scaling exponent of the root mean square displace- 
ment, and D2 the correlation dimension of the multifrac- 
tal eigenstates. We have further provided a conceptually 
new interpretation of D2 as the effective dimension for 
electron motion in the multifractal space. Following this 
interpretation, the above equality can be generalized to d- 
dimensional systems as f3 = /d, as long as the electron 
motion is ballistic in the -space. We have demostrated 
that such an equality applies to exactly solvable 2D and 
3D quasiperiodic systems. When the motion in the " 
space is non-ballistic due to disorder, the equality is re- 
placed by /3 < D^/d, which is in complete agreement 
with all the known results for higher dimensional (d > 1) 
disordered systems. 
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